Solar System experiments do not yet veto modified gravity models 
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The dynamical equivalence between modified and scalar-tensor gravity theories is revisited and it 
is concluded that it breaks down in the limit to general relativity. A gauge-independent analysis of 
cosmological perturbations in both classes of theories lends independent support to this conclusion. 
As a consequence, the PPN formalism of scalar-tensor gravity and Solar System experiments do not 
veto modified gravity, as previously thought. 
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INTRODUCTION 

There are many models in the literature aiming at ex- 
plaining the observed acceleration of the cosmic expan- 
sion discovered with supernovae of type la 0] in con- 
junction with the recent cosmic microwave background 
experiments 2] . One class of models postulates that the 
universe is filled with unclustered dark energy compris- 
ing 70 percent of its energy content. This dark energy 
has exotic properties and, if its effective equation of state 
is truly such that P < —p (where p and P are the dark 
energy density and pressure, respectively), as the obser- 
vations seem to favour [3(|, it is even more exotic and it 
is called phantom energy or superquintessence. Phantom 
energy violates all of the energy conditions and is rather 
difficult to accept because of the possibility of instabili- 
ties, ghosts, and strange thermodynamical behaviour [4j. 
As an alternative to such exotic physics, it has been sug- 
gested that perhaps gravity should be modified at large 
scales 0, by introducing in the gravitational sector 
terms non-linear in the Ricci curvature R. This way, 
one can dispense entirely with exotic dark energy. Apart 
from this ad hoc justification, there are also motivations 
(and corrections) for non-linear gravity from M-thcory 
Scenarios based on this idea are called "modified 
xavity" , "non-linear theories" , or "fourth-order gravity" 
In its simplest form, the action is 
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connection, one obtains a third possibility, metric-affmc 
gravity theories 0, |n| . In the following we consider the 
metric approach to modified gravity but the methods and 
the conclusions apply to the Palatini approach as well. 

We briefly recall the dynamical e quiv alence between 
f(R) gravity and scalar-tensor gravity |l2Lll3U14j (for the 
dynamical equivalence between scalar-tensor theories see 
Ref. By introducing an auxiliary field <j> the action 

||TJ becomes 



S = 



if d f/d(f> 7^ 0. This action integral can be written as 

s = y k / d 4 * v=a mm - v (<i>)] + s (m) , (4) 



where 



m = fr, v^) = <t>^-m. (5) 



This action describes a scalar-tensor theory of gravity 
( [HI E3 — see f° r a review) with Brans-Dicke pa- 
rameter u> = 0. The corresponding field equations are 
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where k = 8nG. The corresponding field equations are 
f'Rab -tg ab = V Q V h /' - g ab □/' + kT^ , (2) 



Trivially, if <fi = R, the action © reduces to QJ. Vice- 
versa, by varying the action (© with respect to <j) an d 
assuming that is independent of <j>, one obtains 



where a prime denotes differentiation with respect to R. 
Whenever f(R) is non-linear in R, the Palatini variation 
treating the metric and the connection as independent 
variables produces field equations that are different from 
@, which arc obtained with the usual Einstein-Hilbcrt 
variation with respect to the metric only. The "Pala- 
tini approach" is widely used in cosmology, in addition 
to the usual "metric formalism" Furthermore, if 

the matter part of the action also depends on the 



(R - 4>) /"(</>) = o , 



(8) 



which yields <p = R provided that /" ^ (a prime 
now denotes differentiation with respect to cf>). Simi- 
larly, one shows that f(R) gravity in the Palatini for- 
malism is equivalent to aw = —3/2 Brans-Dicke theory 
when the matter action is independent of the connection 
(e.g., 0]). This dynamical equivalence between theories 
can be quite useful but it should not be abused. It has 
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been used to constrain f(R) gravity based on Solar Sys- 
tem bounds on the post-Newtonian parameters of scalar- 
tensor gravity [lj, [19, |2fJ . The underlying logic is that 
deviations from general relativity (GR) are not detected 
in our local spacetime neighborhood, therefore these de- 
viations (if they exist) must be small. This "closeness 
of f(R) gravity to GR" implies that f(R) gravity re- 
duces to GR in an appropriate limit, which we address 
here. While there is in principle no problem in taking this 
limit directly in modified gravity, the dynamical equiva- 
lence with scalar-tensor gravity breaks down in this limit. 
In fact, GR corresponds to f(R) = R and /" = 0, while 
the dynamical equivalence requires f" ^ 0. This fact has 
been overlooked and the dynamical equivalence has been 
used beyond its realm of validity in the limit to GR by 
advocating the parametrized post-Newtonian (PPN) for- 
malism of scalar-tensor gravity. This procedure is invalid 
and it explains why opposite claims of compatibility /non- 
compatibility of f(R) gravity with Solar System experi- 
ments occur in the literature - worse, even the Newtonian 
limit is the subject of dispute jl9ll2(l|. 

The limit to GR of the equivalent scalar-tensor theory 
is more general than the weak field limit: it includes the 
strong field regime and it turns out to be a singular limit, 
as shown below. Moreover, the limit of Brans-Dicke the- 
ory to GR in vacuum, or in the presence of "conformal" 
matter (i.e., with vanishing trace of the stress-energy ten- 
sor), is riddled with problems. Therefore, one must be 
particularly careful in basing all of one's conclusions on 
the compatibility with Solar System experiments (which 
test the gravitational field in vacuo or at very low densi- 
ties) solely on the equivalence with scalar-tensor gravity. 
In the following section we discuss the direct limit of f(R) 
gravity to GR without using the equivalence with scalar- 
tensor gravity, and then the corresponding limit to GR 
of the equivalent scalar-tensor gravity, showing the prob- 
lems arising in this last situation. We do not want to 
commit ourselves to specific choices of the function f(R) 
(e.g., the CDTT model 0,0, etc.) but we consider a 
general form of the function f(R). 



THE LIMIT OF f(R) GRAVITY TO GENERAL 
RELATIVITY 



Perhaps the easiest way to consider the limit of f(R) 
gravity to Einstein's theory consists of introducing a 
small parameter e such that f(R) can be expressed as 
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f(R)=R+e<p(R) 



(9) 



The action of GR Son. 
obtained in the limit e 
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0. The field equations 
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which, in the limit e — > 0, formally reduce to the Einstein 
equations G a b — kT^T" . So, there is no problem in taking 
the limit of f(R) gravity to GR directly. Let us consider 
now the "equivalent" scalar-tensor theory J2J. Although 
the conventional way to obtain this limit is letting the 
Brans-Dicke parameter ui go to infinity, here to is fixed to 
be zero. The limit to GR can again be obtained by letting 
e go to zero in the equations of scalar-tensor theory, which 
is equivalent to taking the limit <fi — > constant. Assuming 
that f(R) is given by eq. ©, the non-linear theory JTJ is 
equivalent to © with 



^) = l + e^(0), 



V{cj>) = e{ip'4>-ip) (11) 



provided that /" ^ 0. Now, in the limit e — > 0, 
/" = eip" — > and the equivalence is broken. For- 
mally, the field equation JJjJl reduces to the Einstein equa- 
tion while 10 is identically satisfied. There are however, 
problems with this procedure. One should also find the 
asymptotic behaviour of the field as e — * 0. The situa- 
tion is analogous to the standard limit to GR of Brans- 
Dicke theory, the prototype of scalar-tensor theories, in 
which GR is usually obtained by taking the limit lo — > oo. 
The standard textbook presentation provides the asymp- 
totic behaviour of the Brans-Dicke field 4>bd- 



O 



(12) 



where 0o is a constant [2j|. But when the trace of the 
energy-momentum tensor of matter vanishes, the 

lu — ► oo limit fails to give back GR — this phenomenon 
is reported for a number of exact Brans-Dicke solutions 
pjj and is identified as a general feature of Brans-Dicke 
theory explained by a restricted conformal invariance en- 
joyed by the theory when T< m ) = H3- This anoma- 
lous behaviour is intimately linked with the asymptotics 
displayed by the Brans-Dicke field in these situations 

4>BD = 0O + O \-j=\ (t< W > = 0) . (13) 

Similarly, the examination of the asymptotics of the 
scalar field <j> as the parameter e tends to zero should 
provide a useful check of the limiting procedure. In the 
scalar-tensor equivalent of f(R) gravity, the Brans-Dicke 
parameter is fixed to be to — in the metric formalism 
(lu = —3/2 in the Palatini formalism) and we must nec- 
essarily come up with a different way of taking the limit 
to GR, hence the possibility considered of letting e going 
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to zero while <fi becomes constant. In this case, we should 
obtain a reasonable asymptotic behaviour for the fields 
g ab and 0, say 



9ab = 9 ab +th ab , 



where ffi? is the general-relativistic metric, O is a con- 
stant, and the remainder r(e) tends to zero as e — > 0. 
However, this is not the case. By inserting eq. ijHJ into 
eq. JBJl one obtains 



= O + r(e) , (14) 



G a b — 
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while 
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where the indices are raised and lowered with g3 ■ Fur- 
ther substitution of eq. (|14J) yields, in the limit e — * 0, 



r(e) 



(17) 



The remainder r(e) diverges instead of vanishing: e — > 
is a singular limit of the "equivalent" scalar-tensor ver- 
sion of f(R) gravity while the direct limit e — > of /(i?) 
gravity does not lead to this problem. Again, this re- 
flects the breakdown of the dynamical equivalence in the 
limit to GR in which /" — > 0. Note that the procedure 
employed here parallels the procedure used to obtain the 
estimate (O for Brans-Dicke theory [2! Ijfl . 



STABILITY OF DE SITTER SPACE IN 
MODIFIED AND SCALAR-TENSOR GRAVITY 



In this section we consider the cosmological dynam- 
ics of modified gravity. By assuming the spatially flat 
Friedmann-Lemaitre-Robertson- Walker (FLRW) metric 



ds 2 



-dt 2 + a 2 (t) (dx 2 + dy 2 + dz 2 ) 



(18) 



in comoving coordinates (t, x, y, z), the field equations (0) 
of modified gravity reduce to 



2 _ 1 f f(R)-Rf'(R) ■ 
H - ?777m mRf 



3f'(R) \ 
W m) } , 



(R) 



(19) 



2H + 3H 2 = - 



JJ^Ur) f"(R)+2HRf"(R) 

+Rf"(R) - \ {f(R) - Rf'(R)] + J . 

(20) 



Consider, for simplicity, the situation in which the cur- 
vature terms dominate the dynamics and p( m > and p( m ) 
are negligible. Then, the scale factor a(t) enters the field 
equations only through the Hubble parameter H = a/a 
and it is natural to use H as the dynamical variable. 
The field equations are of fourth order in a (hence the 
name "fourth order gravity" given to f(R) theories) and 
of third order in H . The main result of this section is 
that the u> = scalar-tensor theory does not provide the 
same stability condition derived in f(R) gravity, but two 
different ones according to the type of perturbations con- 
sidered. Therefore, these two theories are inequivalent 
with regard to stability. 

The equilibrium points in the phase space (^H, H, H^j 

are de Sitter spaces characterized by constant Hubble 
parameter Hq given by 



Hi 



A 
6/5 



(21) 



and Rq = 12Hq. The stability of these de Sitter spaces 
with respect to both homogeneous and inhomogeneous 
perturbations was studied in Ref. |3^|, with the result 
that the stability conditions with respect to both types 
of perturbations coincide and, in our notations, are ex- 
pressed by |3J| 



fo fo 



> 



(22) 



The study of stability with respect to inhomogeneous 
perturbations, which are inherently gauge-dependent, re- 
quires the use of a gauge-invariant formalism. It is coun- 
terintuitive that the stability condition with respect to 
inhomogeneous perturbations is not more restrictive than 
the corresponding stability condition with respect to ho- 
mogeneous perturbations. This is not the case, for ex- 
ample, in scalar-tensor gravity. Stability in scalar-tensor 
gravity was also studied in Refs. |U EH- but the cases 
lu = and lu = —3/2 were excluded by the analysis. In 
the following subsections we study the stability of the 
de Sitter equilibrium points in the lu — equivalent of 
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metric f(R) gravity. As expected, the stability condition 
with respect to inhomogeneous perturbations turns out 
to be different from the stability condition with respect 
to homogeneous ones. 



Stability with respect to homogeneous perturbations 
in lo — scalar-tensor gravity 



This is the desired stability condition with respect to 
homogeneous perturbations in the lo = scalar-tensor 
theory that is supposed to be equivalent to f(R) grav- 
ity. This condition is different from IL'l'li showing that, at 
best, the equivalence should be treated with care. The 
stability condition (|33l) should be compared with the sta- 
bility condition with respect to inhomogeneous perturba- 
tions, which we derive in the next subsection. 



In the lo = scalar-tensor theory described by the 
action J3J the field equations become, in the metric l|18|) 
and in the absence of matter, 



1 



3f'H 2 = -(^f-f)-3Hf 



-2f'H = f" [</>) + f"<t> - Hf'q 
f'R - V = . 



(23) 
(24) 
(25) 



Manipulation of eqs. [12: il) and (|24|1 leads to the Klein- 
Gordon-like equation for <p 



<t> + *H<])= — 



-f'R-3f" (V) 2 + 2(0/'-/) 



(26) 

The equilibrium points in this picture correspond to de 
Sitter spaces with constant scalar field given by 



H 2 = 



A 

6/5 ' 



i _ 2 /o _ p 
»o — —p- — -n-o 

JO 



(27) 



Inhomogeneous perturbations in lo — scalar-tensor 
gravity 

The analysis of inhomogeneous perturbations necessar- 
ily requires a gauge-independent formalism. We adopt 
the covariant and g aug e-invariant formalism of Bardeen- 
Ellis-Bruni-Hwang 1421 in the version of Hwang valid for 
generalized gravity |43j. The metric perturbations are 
given by 



5oo 



9oi 



-a 2 (1 + 2 AY) 
-a 2 BY, , 



(34) 
(35) 



g tJ = a 2 [hij (1 + 2H L ) + 2H T Y tj ] , (36) 

where Y, Yi, and ly are the scalar, vector, and tensor 
harmonics, respectively, satisfying 



Homogeneous perturbations of the de Sitter fixed points 
are described by 

H(t) = H + 5H(t) , 4>{t) = 0o + 5<t>{t) , (28) 

and obey the evolution equations 

12iJ /o SH + (6H 2 - <f>o) fH S(b + 6-ff /o 5$ = , (29) 



ViV 4 Y = -k 2 Y 



YiA 



1 



■v<v,-y 



i 



■Yht 



(37) 



(38) 



-2f'8H = fg6$-H f"> 



(30) 



6%5H + 24H f^'8H + (l2H 2 f^ - O / O " - fa) 8</> = 0, 

(31) 

where the constraint <p = R has been used. By eliminat- 
ing SH and using the zero-order equations for de Sitter 
space, one obtains 



k 2 ' " J ~ '3 

Here hij is the three-dimensional metric of the FLRW 
background and the operator Vj is the covariant deriva- 
tive associated with hij, while k is an eigenvalue. The 
gauge-invariant variables used are Bardeen's potentials 
and the Ellis-Bruni variable 



H 



Hr 



Hj 1 



5<j> + 3H 5<j) 



1 



3/a/i'r 



(/o) - 2/o 



. (32) 
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Stability is achieved, and the perturbations 8<p do not run 
away, when the effective mass squared in this harmonic 
oscillator equation is non-negative, i.e., 



Acf> = 6cf) + - , 



B — — Ht 

k 



(39) 



, (40) 



(41) 



(/o) 2 ~ 2/o 



> o 



(33) 



with A/ and AR defined similarly to the last equa- 
tion. The first order equations of motion for the gauge- 
invariant perturbations can be found in Ref. |43| : they 
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simplify considerably in the de Sitter background, reduc- 
ing to 



AV> + 3H Aip - 



^-4ff 2 W + ^Ai? = 0, (42) 



1 ( Aib Aib \ 

-$h + H^ A = - I — Hq , 43 

2 \ ipo ipo I 



fcV k 2 Aib 

a J 2ar ipo 



Qa + ^h = — r~ , 

ipo 



Ht + 3Hq Ht H — 7t Ht — 



Vn 

$ H + 3H Q $ H - H $ A ~—<$>a 

2ip Q 



(44) 



(45) 



(46) 



ipo 



Axp V 



ipo ZiPo 



(47) 



The term (k/a) can be neglected at late times when a 
diverges exponentially fast and the stability condition for 
the perturbations Acp (and = §a oc Acp) is therefore 



Jo Jo 



> 



(52) 



and 



This is the desired stability condition with respect to 
inhomogeneous perturbations. It is different from the 
condition for homogeneous perturbations l|33|l and it co- 
incides with (|22() . Not only one of the stability condi- 
tions with respect to homogeneous and inhomogeneous 
perturbations l|33|) and Q52JI of the ui — scalar-tensor 
"equivalent" of f(R) gravity fails to coincide with the 
condition (|22|) . but they also differ from each other. The 
purported equivalence is not a true equivalence. The rea- 
son should be looked for in the fact that, while in a true 
scalar-tensor theory the Brans-Dicke-like field is a true 
dynamical field whose evolution is only ruled by the dy- 
namical field equations, in the theory considered here <p 
is forced to obey the additional constraint <p = R, thus 
limiting its natural evolution. In other words, we have 
gone from fourth order equations to second order equa- 
tions by adding a scalar degree of freedom to the theory, 
which has now spin content in addition to spin 2, but 
the scalar degree of freedom is somehow constrained by 
the condition Acp = AR. 



DISCUSSION AND CONCLUSIONS 



AR = 6 



2 k 2 

<S> H + <iHo<f> H + -— $ H - H <P A 

3 a z 



^-4HS)<I., 



(48) 



where a = a» e Hot is the unperturbed scale factor. Equa- 
tion i|4T)|) exhibits a positive effective mass squared k 2 /a 2 
for the tensor modes Ht, therefore de Sitter spaces are 
always stable with respect to these modes, to linear order. 
By using eq. Q45JI to eliminate &a and Taylor-expanding 
Aip = ip' A(p + ... , one easily obtains 



H 



$A 



2<Po 



A4> 



(49) 



and 



AR 



A<p + 3H a A<p + — - AH 2 ) Aq) 



(50) 

By using the fact that AR = A(p, the equation for the 
scalar perturbations A<p is obtained: 

A4> + 3H A<b+ (^-AH 2 + ^) AcP = Q. (51) 
\a 2 3% J 



While the limit to GR e — > in the action JIJ of modi- 
fied gravity presents no problems of principle, the limit to 
GR of the "equivalent" scalar-tensor theory is ill-defined. 
In view of the problems presented by this limit in Brans- 
Dicke theory, special care is advised when using the dy- 
namically equivalent scalar-tensor theory to analyse the 
weak field limit of modified gravity. As shown above, the 
limit to GR of the "equivalent" u> = Brans-Dicke theory 
is a singular one, illustrating the fact that the equivalence 
breaks down in this limit. A posteriori it is easy to see 
that this is implied by the fact that f"(R) — * in this 
limit. 

Another issue is the following: if the dynamical equiv- 
alence were to hold in the limit to GR, the experimental 
bound to > 40000 01 would be in violent conflict with 
the values of the Brans-Dicke parameter u> — or —3/2 
obtained, unless the field <p is short-ranged. The effective 
mass of <f> is given by 



l eff 



= V" = e(v" + <Pv"') 



(53) 



and vanishes as e — > 0, making it impossible to suppress 
the violation of the bounds on cj, for any form of the 
function f(R). This contradicts the results of Refs. 20, 
45] which support the viability of the weak field limit of 
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the theory for specific forms of f(R). As it turns out, the 
PPN limit of the associated scalar-tensor theory bears no 
relation to the weak field limit of the original modified 
gravity theory. 

As a consequence, the conclusion 0, ^| that mod- 
ified gravity always violates the stringent Solar System 
bounds on scalar-tensor gravity is invalid. The is- 
sue of the correct Newtonian and post-Newtonian limit 
of such theories is still open and should be approached 
directly without using the dynamical equivalence discov- 
ered in Refs. 0, El , which is still useful for other pur- 
poses. The regime that is more interesting, however, is 
the one in which the deviations from GR in the Solar Sys- 
tem are small. A complete study of the PPN limit of gen- 
eral modified gravity scenarios without resortin g to the 
equivalent scalar-tensor theory (initiated in Refs. [2(1 kg ) 
will be presented elsewhere. 

Another apparent problem with the limit to GR lies in 
the fact that (f> must become constant: because cf> = R, 
were this limit correct, it could only reproduce solutions 
with constant Ricci curvature (which includes vacuum 
solutions and solutions sourced by conformal matter). 
Although this could work for vacuum solutions used to 
describe Solar System experiments, it is by no means 
acceptable to have a limit to GR valid only for special 
solutions: the limit must apply to the general theory. 
However, we believe that this second problem is not very 
relevant because, when /" vanishes in eq. (jHJ) in this limit, 
the equality <fi = R is no longer enforced. 

Finally, the viability of modified gravity scenarios does 
not rely only on its correct weak field limit: other issues 
are the presence of ghosts and instabilities and, of course, 
a correct cosmological dynamics. These have been con- 
sidered separately in the literature 
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